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ABSTRACT 

Recent work using virial mass estimates and the quasar mass-luminosity plane has 
yielded several new puzzles regarding quasar accretion, including a sub-Eddington 
boundary on most quasar accretion, near-independence of the accretion rate from 
properties of the host galaxy, and a cosmic synchronization of accretion among black 
holes of a common mass. We consider how these puzzles might change if virial mass 
estimation turns out to have a systematic bias. As examples, we consider two recent 
claims of mass-dependent biases in Mgll masses. Under any such correction, the sur- 
prising cosmic synchronization of quasar accretion rates and independence from the 
host galaxy remain. The slope and location of the sub-Eddington boundary are very 
sensitive to biases in virial mass estimation, and various mass calibrations appear to 
favor different possible physical explanations for feedback between the central black 
hole and its environment. The alternative mass estimators considered do not simply 
remove puzzling quasar behavior, but rather replace it with new puzzles that may be 
more difficult to solve than those using current virial mass estimators and the Shen 
et al. (2008) catalog. 

Key words: black hole physics — galaxies: evolution — galaxies: nuclei — quasars: 
general — accretion, accretion discs 
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1 INTRODUCTION 

Recent work developing the quasar mass-l uminosity plane 



from a combination of large quasar catalogs JSchneider et al 



201C ) and virial mass estima tors l|Vestergaard fc Peterson 



!( 1 ^^^^^^^^^^^^^^^^^^^^^^^^^^^ 

McLure fc D unlop 2004) has provided a new picture of 



Type 1 quasar accretion seemingly i ncons istent with existing 
models jsteinhardt fc Elvisll2010bl . l201ll ). including: 

1. Characteristic Accretion Rate at Fixed M and 

z: At any fixed combination of mass and redshift. Type 1 
quasar accretion is apparently only possible within a narrow 
(< 1 dex) luminosity range around a central, characteristic, 
sub-Eddington accretion rate. Therefore, ephemeral proper- 
ties of host galaxies such as star formation rate and recent 
merger activity appear to have a minimal effect upon the 
accretion rate of the central supermassive black hole. This 
appears to be in sharp contrast with Seyfert galaxies and 
other active galactic nuclei with low Eddington ratios, for 
which properties of the host galaxy are more important. 

2. Accretion is Time-Dependent: The characteristic ac- 
cretion rate is time-dependent at fixed mass, decreasing to- 
wards lower redshift. Even though galaxies of the same stel- 
lar mass may virialize at different times and have different 
morphologies at the same redshift, the accretion rates of 
their central supermassive black holes are well synchronized. 



3. Sub-Eddington Accretion: The characteristic accre- 
tion rate is near the Eddington luminosity for the lowest- 
mass quasars at every redshift, but at higher masses falls 
increasingly short of Eddington. Quasars at most combi- 
nations of mass and redshift are constrained to lie below 
Eddington, with the highest mass quasars (and most of the 
mass growth) taking place closer to L/LEdd ^ 0.1. 

There appears to be a universality to quasar accretion 
implying that only a few, simple physical processes may play 
important roles. If correct, these results may be incompati- 
ble with our current understanding not just of quasar accre- 
tion and turnoff but perhaps also of structure formation. 

Thes e analyses were pe rformed using virial masses pro- 
vided by ih en et al.l (|2008l ) in an extension to the Sloan 
Digital Sky Survey (SDSS) quasar catalog (Schneider 2007i). 
Since the publication of these masses, there have been sev- 
eral investigations into whether virial masses might sys- 
tematically mis-measure black hole masses, particularly for 
Mgll-based masses. Perhaps the existing virial mass scal- 
ing relationships must be corrected to account for quasar 
wind, thermal effects, or other non-virial motion in the 
quasar broad-line region (BLR) lOnken fc Kollmeieij l2008l : 
iRisaliti et al] l2009l : iMarconi et all |2009| ). Alternatively, at 
least the statistical uncertainty in virial mass estimation ap- 
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pears to be better than previously believed, and might even 
be better t han that of more difficult r e verberation mapping 
techniques (IStei nhardt fc Elvisll2010d ). lRafiee fc Halll (|20ld i 
claim that different methods of determining the velocity of 
gas in the quasar BLR show evidence of a mass-dependent, 
systematic error. They also claim that correcting for this 
error returns all quasars to lying at their Eddington lumi- 
nosity. 

These reports of possible biases in virial mass es- 
timation must also lead us to re-examine the con- 
clusions drawn from the quasar mass- luminosity plane. 
Virial mass estimation is inherently a difficult proposi- 
tion, and estimators are calibrated against only 35 low- 
redshift q uasars with reverberation mapping-based mass 
estimates ([Vestergaard fc Petersonll2006l : Iwang et al]|2009l : 



iPeterson et all l2004l l. It is therefore essential to under- 
stand how s ensiti ve each of the conclusions drawn using the 
ih en et all ||2008D mass catalog are to biases in virial mass 
estimation. 

In general, potential biases can be divided into four cat- 
egories: 

(i) Afflne Biases: Biases that affinely translate virial 
mass estimate, shifting, stretching, or contracting the dis- 
tribution in the log M — log L plane. While this may seem a 
restrictive definition, because the velocity of broad-line gas 
is not strongly correlated with its radius, most currently- 
proposed corrections to virial mass relationships are affine. 

(ii) Correctable Biases: Biases that do not affect virial 
mass estimates uniformly across a large catalog, but for 
which each spectrum contains enough information to cor- 
rect them. For example, a hypothetical bias might only af- 
fect quasars with strong Hell lines, but could be corrected 
by properly subtracting Hell before fitting the H/3-OIII com- 
plex. 

(iii) Uncorrectable Biases: Biases that affect different 
objects in different ways and cannot be corrected from the 
individual spectra, but can be understood well enough to 
be modeled in a probabilistic way for large catalogs. For 
example, a hypothetical bias might arise from a time de- 
lay between variability in the quasar continuum and re- 
sponse of the broad-line region, leading to errors when using 
monochromatic continuum luminosity to estimate the BLR 
radius in some objects. 

(iv) Byzantine Biases: Along the lines of IPease et al.l 
arbitrary biases maliciously chosen as a worst-case 

scenario; one might use these to determine what we have 
learned with absolute certainty. 

In this work, we will consider the effects of affine bi- 
&SGS ^ cL catego ry that includes mo st cu rrently-proposed cor - 
rections. Both lRisaliti et all (120091 ') and lRafiee fc Halll (|2010l ') 
propose systematic, linear corrections for Mgll-based es- 
timators relating the true black hole mass Mbh and the 
ShenOS mass Mshen as Mbh = AMshen -\- B. A compres- 
sive correction with A > 1 will act to narrow the quasar 
meiss distribution from its appearance in the ShenOS cata- 
log, while an expansive correction with ^4 < 1 will broaden it. 
It is instructive to consider the effects of both types of cor- 
rections, biasing the ShenOS masses in different directions, 
on the key puzzles emanating from the M — L plane analysis 
of the SPSS catalog. In § p we c onsider the proposed cor- 
11( 20101 ). a compressive correction 
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rection from lRafiee fc Hal 



203), 



with yl < 1. In §|3] the implications of an e xpansive correc- 
tion with A > 1, from iRisaliti et al] l|2009l ). are discussed. 
Finally, § |3] considers the lessons learned from this analysis. 



2 A COMPRESSIVE CORRECTION (RAFIEE 
& HALL 2010) 

iRafiee fc Halil (|2010l ) (RHIO) consider the possibility that 
the the full-width half-maximum (FWHM) of quasar broad 
emission lines might not directly indicate the velocity of 
gas in the broad-line region, i.e., that motion in the broad 
line region may not be predominantly virial. They instead 
fit Mbh oc FWHMif^'*° using the MLIN-IX_ERR fit- 
ting method of lKellvl (|2007l ). Although this possibility might 
mean aban doning the virial assumptio n used in virial mass 
estimation (|Steinhardt fc Elvisll2010al ). we use it as an ex- 
ample of a possible mass-dependent bias that produces a 
compressive correction, narrowing the mass distribution at 
fixed redshift. 

Applying the Rafiee fc Hall (RH) formula produces ap- 
proximately log Mrh = 0.67 log Mshe„ -f 3.0. Although the 
RHIO correction itself is only calibrated for Mgll, since virial 
mass estimators using different broad emission lines are cal- 
ibrated against each other, we will consider the effects of 
this adjustment at all redshifts, not just 0.4 < z < 2.0 where 
Mgll is available. A compressive correction such as this does 
not remove all of the puzzling behavior in the quasar mass- 
luminosity plane, but does recast some puzzles in a different 
light. 

2.1 Puzzling Behavior 

Specifically, this adjustment produces the mass-luminosity 
plane shown in Figure [T] in 12 redshift bins ranging from 

0. 2 < z < 4.1. These distributions exhibit the following 
substantial puzzles: 

1. Characteristic Accretion Rate at Fixed M,z: At 
any fixed combination of mass and redshift, the quasar lu- 
minosity range spans a range of less than O.S — 1 dex. Because 
SDSS views objects across nearly half the sky, the catalog 
includes quasars that should lie in host galaxies with very 
different ephemeral properties such as star formation rate 
and recent merger activity. Thus, these ephemeral proper- 
ties can only be important at the < 1 dex level to the quasar 
accretion rate. This puzzle also appears when using ShenOS 
masses. 

2. Linear Mass Dependence: At fixed redshift, the char- 
acteristic accretion rate increases nearly linearly with mass. 
In the ShenOS picture, the characteristic accretion rate in- 
creased between L oc Afg^^""'^, and thus for log Mrh ~ 
0.67 log Afshen -1-3.0, the characteristic accretion rate will be 
closer to L (X M. This is best demonstrated by the behav- 
ior of the sub-Eddington boundary, or the luminosity of the 
most luminous quasars at each mass, using different mass es- 
timation techniques (Fig. [2| . The sub-Eddington boundary 
should run parallel to the characteristic accretion rate, but 
is easier to measure precisely because it is determined us- 
ing the brightest and best-measured quasar spectra at each 
mass and redshift. Because the slope appears to vary with 
redshift, the accretion rate cannot be exactly linear in M 
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Figure 1. A contour plot of the mass-luminosity distribution in 12 different redshift ranges (0.2 - 4.1) using two mass estimators (H/3 and 
Mgll) using the Rafiee & Hall (2010) mass estimators. Quasars at decreasing redshift lie increasingly below their Eddington luminosity 
(black, dashed). 



(or at any other constant value) at all redshifts under any 
time-independent adjustment to the ShenOS masses. Never- 
theless, the accretion rate remains nearly linear. Since the 
Eddington luminosity is also Lsdd oc M, this means quasars 
at fixed redshift all lie at nearly the same Eddington ratio 
L/LEdd, rather than at different L/LEdd oc Mgl^^'^~^'^\ 

3. Accretion is Time-Dependent: The characteristic ac- 
cretion rate at fixed mass, or chararacteristic Eddington ra- 
tio across all masses in the RIflO picture, increases with 
increasing redshift (Fig. [!}. Near z ^ 2, quasars are typ- 
ically able to reach but not exceed their Eddington lumi- 
nosity. However, at lower redshift, quasars all miss Edding- 
ton by increasingly large margins, and at 2 > 3.0, CIV- 
based masses correspond to super-Eddington accretion. Al- 
though CIV masses have a substantially higher statistical 
uncertainty than Mgll (jSteinhardt fc Elv is 2010c), suffer 
from quasar wind and ra diation pressure-induced broaden- 
ing (|M 

arconi et al.ll2009l ). an d have a lar ge scatter when 
compared against Mgll masses l|Shen et al.| :2008). there is no 
systematic disagreement between Mgll and CIV masses in 



the ShenOS sample. Therefore, corrections inducing a mass- 
dependent shift in ShenOS Mgll masses should also apply to 
CIV masses, and the resulting statistically significant super- 
Eddington accretion is either evidence against such a mass 
correction or yet another puzzle. 

4. Narrow Active Mass Range: At fixed redshift, the 
observed range in quasar luminosities spans ~ 1 — 1.5 dex, 
while ShenOS masses span 2 — 2.5 dex. A correction narrow- 
ing the mass range to the point that quasars all lie near a 
fixed Eddington ratio must also result in a mass range span- 
ning 1 — 1.5 dex at fixed redshift. There are three major 
contributions to the width of this mass distribution: 

• Mass Growth: From when they first enter the SDSS 
catalog, quasars grow in mass via luminous accretion. 
Quasars must enter the catalog at least as the lowest-mass 
objects at the redshift they are observed, and must cease 
that luminous accretion at a mass no higher than the max- 
imum mass observed at some later redshift. Since the maxi- 
mum mass is decreasing with redshift, the turnoff mass for 
a quasar in the SDSS catalog must actually be lower than 
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Figure 2. Slope of the sub-Eddington boundary as a function 
of mass for masses estimated using four virial mass calibrators: 
H/3 (green) and Mgll (blue) masses from Shen et al. (2008) are 
compared with Mgll masses from Rafiee &; Hall (2010) (red) and 
Risaliti, Young, & Elvis (2009) (magenta). Only the Rafiee & Hall 
masses yield L oc M (black, dashed) at any redshift, although no 
mass estimate is consistent with this dependency at 2 ~ 1.9. 

the maxim um mass at the redshift it is observ ed. The Soltan 
argument l|Solta nl ll982l : IYu fc Trema"in3l2002l ) demonstrates 
that most supermassive black hole mass is accreted as a 
quasar; the mass growth in this phase must be at least ~ 1 
dex. 

• Different Timing: If all quasars first appeared as the 
lowest-mass quasars at some redshift, the width of the mass 
distribution might correspond to the amount of mass growth 
in a quasar phase. However, if quasars first appear at higher 
masses, they will still lie in the observed mass distribution 
as long as they are quasars, and thus the mass distribu- 
tion will be wider than the growth of any individual object. 
The evolution of quasar host galaxies of a common halo 
or st ellar mass is not well synchronized (cf. (jSpringel et al.l 
|2005| )). yet they will eventually produc e post-turnoff super- 
mass i ve black holes of the same mass (|Ferrarese fc MerrittI 
I2OOOI : iGebhardt et~aLll2000l ) . The lowest-mass quasars in the 
SDSS catalog decrease with decreasing redshift due to de- 
tector sensitivity. Thus, if quasars in these host galaxies all 
turn on at a common mass but at different times, some must 
not be the lowest-mass quasars when they enter the SDSS 
catalog. Differences in the timing of quasar phases will thus 
lead to a mass distribution broader than one due to mass 
growth alone. 

• Mass Errors: Uncertainties in quasar masses will 
artificially broaden the quasar distribution; the observed 
mass distribution is the physical mass distribution con- 
volved with the mass error function. Statistical errors 
alone in virial masses have b een estimated at ~ 0.4 dex 
l|Vestergaard fc Petersonll2006l ). 

For the Soltan argument to apply, luminous accretion 
must comprise the last > 1 dex of mass growth. Thus, the 
timing of quasar phases in different host galaxies of the same 
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Figure 3. The quasar mass-luminosity distribution for all quasars 
at 0.4 < 2 < 2.0, using Mgll masses as described by Rafiee & Hall 
(2010). Quasars at log M/A/q < 8.5 lie below Eddington, with the 
shortfall increasing towards lower mass. 

halo mass must be very well synchronized and mass uncer- 
tainty must have substantially smaller than an ~ 0.4 dex 
scatter in order to produce a mass distribution spanning 
1 — 1.5 dex. Each of these is a puzzle in its own right. 

5. Some Masses Never Reach Eddington: In the 

ShenOS catalog, with IVestergaard fc PetersonI (|2006l ) and 
iMcLure fc Dunlod l|2004l 'l mass estimators, at every redshift 
there is some mass which quasars reach but do not exceed 
their Eddington luminosity, and at every mass there is some 
redshift at which quasars lie at Eddington. With the Rafiee 
& Hall (2010) compressive correction, quasars below z ~ 1.7 
do not reach Eddington. Since only a narrow mass range is 
active at each redshift, quasars with lower masses are never 
seen at Eddington (Fig. [3]). If, as we are assuming in this 
section, this is not a calibration problem, then it is a new 
puzzle requiring a physical explanation. 

2.2 Physical Interpretation 

Like previou s mass estima t ors, t he compressive correction 
proposed by iRafiee fc HaUl tOldt ) produces a quasar mass- 
luminosity plane exhibiting several puzzling behaviors. How- 
ever, some of these features are new to this calibration, 
particularly the near-linear dependence of the characteris- 
tic accretion rate on black hole mass. It is natural to believe 
that the characteristic accretion rate is due to feedback be- 
tween the central black hole and its environment. A linear 
dependence of accretion rate on black hole mass is typical 
of a model in which gravity is dominant. Further, x-ray bi- 
naries, another type of accretion disk system, are capable 
of radiatively-efficient accretion only in a Eddington ratio 
range spanning ~ 1 dex; accretion is radiatively inefficient 
both above and below that range. Using these mass estima- 
tors, quasars, too, can only accrete in an Eddington ratio 
range spanning ~ 1 dex, although the Soltan argument pre- 
cludes radiatively inefficient accretion at a higher Eddington 
ratio being common. 

If the mass-dependence appears easy to understand, 
though, the redshift-dependence is quite difficult. If the ac- 
cretion rate is set purely by gravity, and feedback from the 
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host galaxy is not important, novel high-energy astrophysics 
or physics will be needed to explain this time-dependence. 

Another puzzling feature using ShenOS masses is the 
synchronization of quasar accretion in host galaxies of the 
same halo mass. A compressive correction appears to require 
even stronger, near-perfect synchronization in order to re- 
main consistent with the Soltan argument. This picture re- 
quires either very well synchronization of the first epoch of 
star formation in each of these host galaxies or, more likely, 
some sort of primordial black hole seeding mechanism. 

3 AN EXPANSIVE CORRECTION (RISALITI, 
YOUNG, & ELVIS 2009) 

iRisaliti et al] (|2009l ') (RYE09) calibrate Mgll-based masses 
against H/3, finding that the Mgll masses systematically dis- 
agree with H/3 masses as lo g Mgg = 1.8 log MmqII — 6.8. As 
discussed in [Stcinhardt & Elvij (|2010d ). this may be due to 
difficulties in separating the Mgll line and continuum from 
nearby Fe lines, and as such, the correction may be sub- 
stantially smaller or unnecessary for the brightest quasars 
with the highest signal-to-noise spectra at each redshift. Evi- 
dence for this includes good agreement in the sub-Eddington 
boundary slope calculated using unadjusted and Mgll 
masses at a common redshift, but poor agreement between 
and adjusted Mgll masses (Fig. [2]). The best agreement 
in sub-Eddington boundary slopes comes from a substan- 
tially sm aller, yet also expansive, co rrection to Mgll mass es- 
timates l|Steinhardt fc ElvisllioiM ). We will use the RYE09 
calibration as an example of a possible mass-dependent bias 
that produces an expansive correction, widening the mass 
distribution at fixed redshift. 

3.1 Puzzling Behavior 

Specifically, this adjustment produces the mass-luminosity 
plane shown in Figure |3] in 12 redshift bins ranging from 

0. 2 < z < 4.1. As with the compressive correction considered 
in § (2] this does not remove all of the puzzling behavior in 
the quasar mass-luminosity plane, but rather includes the 
following surprises: 

1. Characteristic Accretion Rate at Fixed M,z: At 
any fixed combination of mass and redshift, the quasar lumi- 
nosity range spans a range of less than 0.8—1 dex. As before, 
since SDSS quasars should lie in host galaxies with very dif- 
ferent ephemeral properties such as star formation rate and 
recent merger activity. Thus, these ephemeral properties can 
only be important at the < 1 dex level to the quasar accre- 
tion rate. This puzzle also appears when using Shen08 and 
RHIO masses, and can only be removed by a correction that 
interprets quasars with different Shen08 masses as having 
identical black hole masses. 

2. Sub-linear Mass Dependence: At fixed redshift, 
the characteristic luminosity increases with mass as L cx 
^ry~e' '^' with the slope closer to 0.4 at most redshifts (Fig. 
[2)). This dependence is very sensitive to mass-dependent 
biases, pushing towards a linear mass dependence for ex- 
tremely compressive corrections and towards only slight 
mass dependence for extremely expansive corrections. 

3. Accretion is Time-Dependent: The characteristic 



accretion rate at fixed mass, or chararacteristic Edding- 
ton ratio across all masses in the RYE09 picture, increases 
with increasing redshift (Fig. For example, quasars with 
log Mrye/Mq = 10**'^ reach Eddington at z ~ 2, but lie in- 
creasingly sub-Eddington at lower redshifts. Even the most 
luminous quasars at log Mrye/Mq = lO^ '^ and z ~ 0.2 lie 
nearly 1 dex below their Eddington luminosities. As with 
the existence of a characteristic accretion rate, this puzzle 
also appears using Shen08 and RHIO masses; it can only be 
removed with a time-dependent mass adjustment. This may 
be the least surprising puzzle, as it could be described as 
quasar downs i zing a nalogous to galactic and star formation 
iMadau et al.l l|l998f ) downsizing. 

4. Extraordinarily Massive Black Holes: The most 
massive known supermassiv e black hole, OJ 287, has a mass 
estimated at 1.8 x Mq l|Valtonen et al.ll200^ . Applying 
the RYE09 adjustment to Mgll masses produces a more mas- 
sive population at 2 > 1.6 (including log Mbh/Mq > 10.5), 
with the number of more massive objects increasing to- 
wards higher redshift. Although some of these masses use 
higher-uncertainty CiV-based estimates, at 1.6 < z < 2., 
these extraordinarily massive black holes ar e found using 
Mgll masses, which are directly calibrated bv IRisaliti et al.l 
(2009). Such immense black holes would need to lie in im- 
mense galaxies in order to lie on the M — a relation today. 

5. No Eddington Accretion at z ~ 1: At 2; ~ 2 and 
z < 0.8, the least massive quasars reach their Eddington lu- 
minosity, while more massive quasars fall increasingly short. 
This is an a djusted version of the sub- Eddington boundary 
described in ISteinhardt fc ElvisI ([2010b). However, between 
these redshifts, all quasars, including those with the low- 
est black hole masses, fall short of Eddington. Because the 
lowest-mass quasars also lie at luminosities near the SDSS 
detection limit, this may merely be a selection effect. 

3.2 Physical Interpretation 

Like other mass estimators, th e expansive correction pro- 
posed by IRisaliti et al. I l|2009t ) produces a quasar mass- 
luminosity plane exhibiting several puzzling behaviors. Some 
of these behaviors, including the narrow luminosity range at 
fixed mass and redshift and the time-dependence of that 
characteristic accretion rate, are common to all virial mass 
estimators. A unique feature of this correction is an accre- 
tion rate L cx 

The accretion rate is likely due to interaction between 
the central black hole and its environment, and a sub-linear 
mass-dependence likely indicates a balance between gravi- 
tational and non-gravitational effects. One model with ap- 
proximately the right mass dependence considers a compe- 
titio n for infalli ng gas between accretion and star formation 
(Th ompson et a l. 2005), yielding L cx M^/^. 

As before, although models may exist for the mass- 
dependence, the time-dependence is more difficult. The laws 
of general relativity are time-independent, while a black hole 
has no memory; just mass, charge, and spin. We might ex- 
pect the feedback mechanism to be sensitive to properties of 
the host galaxy, but because accretion is well-synchronized 
between host galaxies containing supermassive black holes of 
a common mass, ephemeral properties of the galaxy cannot 
be critical. As in the Shen08 and RHIO pictures of quasar 
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Figure 4. A contour plot of the mass-luminosity distribution in 12 different redshift ranges (0.2 < z < 4.1) using the Risaliti, Young, & 
Elvis mass estimators for Mgll masses and a corresponding adjustment for H/9 and CIV masses. Higher-mass quasars at fixed redshift 
lie increasingly below their Eddington luminosity (black, dashed). 



accretion, the time-dependence continues to defy explana- 
tion. 

4 DISCUSSION 

We have considered the effects of affine biases in virial mass 
estimation by exa mining two recent r e ports of biases in di f- 
ferent directions (|Risaliti et al.ll2009l : Rafiee fc Halll |2010| ). 
Correcting for one would expand the mass distribution at 
fixed redshift, while the other would compress it. Common 
to all of these virial mass estimators are two major puzzles: 
(1) that quasars at a common black hole mass and redshift 
have the same, characteristic luminosity and (2) that the 
characteristic luminosity is time-dependent. 

It should be immediately apparent that a purely mass- 
dependent error cannot be responsible for these two results. 



Any linear (or more complex) correction to black hole masses 
in which objects at a common ShenOS mass have a common 
black hole mass while objects at different ShenOS masses 
have different black hole masses will yield this same pair of 
puzzles. A related puzzle is that using ShenOS masses, the 
lowest-mass quasars at each redshift are found near their 
Eddington luminosities, while the highest- mass quasars all 
fall short of Eddington. iRafiee fc H all (2010) claim to have 
removed this "sub-Eddington boundary" (SEB) with their 
systematic correction, but because the characteristic lumi- 
nosity is time-dependent, it is, again, impossible for a purely 
mass-dependent correction to Mgll masses to remove this 
boundary at all choices of redshift. 

Another puzzling feature of the quasar M — L distri- 
bution is the mass-dependence o f the characteristi c lumi- 
nosity at fixed redshift. Using the IShen et al.l (|2008l ) exten- 
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sion of t he SPSS cata log based upon previous mass cali- 
brations jMcLure fc Jar vis 200^; iMcLure fc Dunlod |2004 
IVestergaard fc PetersonI 12009 ) produces Lc cx M°-*-°-'' , 
which is presently unmatched by physical modeling. Unlike 
the previous two puzzles, the slope of the ma ss-dependence 
is hig hly sensitive to mass correction. The iRafiee fc Halll 
(|20ld ') correction produces a nearly linear dependence, im- 
plying that the characteristic accretion rate might be ex- 
plained due purely to gravity from the central black hole. 
Such an explanation would help us to understand why the 
characteristic accretion rate is universal, and does not de- 
pend upon properties of the host galaxy such as star forma- 
tion rate, me r ger a ctivity, and age since virialization. The 
iRisaliti et al.l (|2009! ) correction, on the other hand, is a good 
match for a feedback model describing competition between 
gravitational attraction and star formation for infalling gas, 
which again might be independent of ephemeral properties 
of the host. 

However, neither mechanism can produce all of the the 
peculiar boundar ies in the M — L plane obtained using 
IShen et all l|2008l ) masses, most importantly in regard to 
the puzzling synchronization of quasar accretion. Moreover, 
both of these mass corrections introduce new, puzzling fea- 
tures in the M — L plane. Using the Rafiee fc Hall correction, 
all quasars will fall short of Eddington at some redshifts, all 
quasars at some mass will never reach Eddington, and the 
quasar distribution may be so narrow that the Soltan argu- 
ment requires perfect synchronization and tiny mass errors. 

Using the Risaliti, Young, & Elvis correction, again all 
quasars fall short of Eddington at some redshifts, and many 
quasars contain black holes more massive than any cur- 
rently known. The Shen et al. masses, on the other hand, 
provide some quasars reaching but not exceeding Edding- 
ton at all masses, provide some quasars reaching but not 
exceeding Eddington at all redshifts, and produce masses 
lying within the bounds of observed post-turnoff super- 
massive black holes near redshift zero. The uncorrected 
masses appear to be the best of the three options consid- 
ered, although a smaller c orrection such as that found in 
ISteinhardt fc Elvij (|2010cl ) by comparing the H/3 and Mgll 
views of sub-Eddington boundaries may also be plausible. 

4.1 Further Remarks 



One of the reasons for skepticism about virial masses has 
been that if taken at face value, the interpretation of 
the mass-luminosity plane appears inconsistent with what 
would seem to be obvious properties of supermassive black 
hole growth. In particular, supermassive black holes are 
composed of gas and dust from their host galaxy, live in 
their host galaxy , and the mass-stellar ve l ocity dispersion 
(M- a) relation (|Ferrarese fc Merrittll2000l : lGebhardt et al.1 
l200d) and black hole mass-galactic luminosity relation 
( Ferrarese fc Fordll2005l ) require that SMBH cease their ac- 
cretion at a final mass dynamically linked to their host 
galaxy. It also seems natural that episodes of stronger ac- 
cretion might be linked to galactic processes such as star 
formation and mergers that might drive gas and dust into 
their central regions. 

It should be noted that although these puzzles re- 
garding SMBH accretion are new, several other evolution- 
ary phases remain poorly understood. There are a wide 



variety of proposed quasar seeding mechanisms jVolonteril 
none of which has yet been supported by direct 
observational evidence. The bla ck hole ma ss-stcUar vcloc- 
ity dispersion {M — a ) relation (|Ferrarese fc Mcrritt 200^; 
iGebhardt eTail I2OO0I) and black hole mass-galactic lumi- 
nosity relation (|Ferrarese fc Fordll2005l ) require that SMBH 
cease their accretion at a final mass dynamically linked to 
their host galaxy, but this link is not understood theoreti- 
cally. In short, including these new results, we must conclude 
that we do not understand how SMBH are born, we do not 
understand how they grow, and we do not understand how 
they die. 

An inability to understand observational results within 
the framework of a physical model should indeed be wor- 
rying, but in this case, the alternative explanations are no 
better. The Sloan Digital Sky Survey quasar catalog ap- 
pears complete when compared aga inst othe r quasa r cat- 
alogs (|Jester et all 120051 : iSteinhard't fc Elvid l2010bl ). and 
while the interpretation of virial masses as indicative of the 
central black hole mass may be in question, our ability to 
measure the component continuum flux and emission line 
FWHM is not. Thus, in order to explain these observations, 
one of two things must be true: 

• Virial masses are not indicative of the central 
black hole mass. Therefore, to explain the observed quasar 
distribution in M, L, and z, a physical model for accretion is 
needed in which the BLR velocity is very tightly correlated 
with the bolometric luminosity, but that correlation is also 
sharply redshift-dependent. 

• Virial masses do measure the central black hole 
mass. Therefore, to explain the observed quasar distribu- 
tion in M, L, and z, a physical model is needed in which 
quasars at fixed mass and redshift truly lie in a narrow 
range in luminosity, and that narrow range is both mass- 
and redshift-dependent. 

It is unclear which interpretation, if correct, would be 
more intriguing. If the BLR velocity is a sufficiently good in- 
dicator of bolometric luminosity, a predictive physical model 
might allow the use of quasars as a bright, high-redshift stan- 
dard candle. If virial masses indeed measure central black 
hole masses, the observed synchroniation of quasar evolu- 
tion continues until turnoff, at which point the SMBH must 
lie on the M — a and M — L relations. Thus, quasars ly- 
ing in two galaxies of the same mass and same redshift, but 
with different star formation rates, different morphologies, 
different merger histories, and different virialization times 
have the same accretion rate. So, we are forced to choose 
between two difficult model-making problems, both requir- 
ing a new understanding of quasar accretion, and neither 
accommodated by existing models. 
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